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Abstract 

We calculate the negative integer moments of the (regularized) characteristic polynomials 
of N x N random matrices taken from the Gaussian Orthogonal Ensemble (GOE) in the 
limit as TV — > oo. The results agree nontrivially with a recent conjecture of Berry & Keating 
motivated by techniques developed in the theory of singularity-dominated strong fluctuations. 
This is the first example where nontrivial predictions obtained using these techniques have 
been proved. 

1 Introduction 

Let H = H T (we here use the symbol T to denote matrix or vector transposition and * to denote 
complex conjugation) be an Nx N random symmetric matrix with real entries distributed according 
to the standard joint probability density of the Gaussian Orthogonal Ensemble (GOE) of random 
matrix theory - 

V{H) = C N e-^\ (1) 

with respect to the measure dH = YliLi dHu Yli<j dHij, where the normalization constant Cn is 
given by 

jV(Af + l) 

1 / N \ 1 

Cn= 2^{^p) (2) 

Z N (n) = dsb(jilN-H) (3) 

denote its characteristic polynomial. We shall here be interested in the negative integer moments 
of \Z\, defined by averaging over the GOE, when Im/i > 0, in the limit as N — > oo. (The positive 
moments of the characteristic polynomials of random unitary-symmetric matrices were calculated 
in pH ; for the positive integer moments it was confirmed in || t hat, as expected, these results also 
apply to the large N limit of matrices in the GOE; see also ]16[.) 

Berry & Keating || (hereinafter referred to as BK) have recently put forward a general conjec- 
ture about the asymptotics of the negative moments of the characteristic polynomials of random 
matrices in the limit as the matrix size tends to infinity when Im^t is scaled by the mean eigenvalue 
density and tends to zero. This conjecture applies to all negative moments, rather than just to 
negative integer moments, and covers all three of the classical random matrix ensembles (i.e. the 
unitary, orthogonal and symplectic ensembles). It predicts a highly non-trivial dependence of the 
asymptotics on the power to which the polynomial is raised. This is in contrast to the case when 



and let 



1 



the large-matrix limit is taken without scaling Im/i by the mean level spacing; then the moment 



In the case of the Gaussian Unitary Ensemble (GUE) of random matrices, the conjecture given 
in BK agrees with the values of the negative integer moments calculated by Fyodorov in || and 
shown to be universal (in the sense that they apply to all unitary-invariant ensembles of Hermitian 
matrices) in However, these values also happen to coincide with the corresponding ones 

when Imp, isn't scaled, so this cannot be said to constitute a test of the non-trivial aspects of the 
conjecture. 

For the GOE of random matrices the conjecture in BK is that the ensemble average of |Z/v(/i)|~ fc 
diverges like e~ vtyk \ as e, Im/x scaled by the mean eigenvalue density, tends to zero, with 



It was suggested in BK (page L4) that, in the notation of the present paper, when k is an integer 
"it is possible that the leading-order power-law behaviour ([|) is multiplied by a power of log . 

Our first aim here is to extend the heuristic arguments developed in BK to recover the log- 
arithmic factor when k takes integral values; this turns out to be simply log ~ for each k. Our 
second aim is then to prove the resulting expression by a direct evaluation of the GOE average. 
In fact, we are able to go significantly further in that we calculate the precise asymptotic form of 
the moments in the appropriate limit. The general expression we obtain (see (^0|) and (|4l|)) takes 
the form of a multiple integral and is interesting in its own right, in particular in view of recent 
endeavours to understand the analytic structure behind the so-called replica limit k —* [T^] . 

The heuristic arguments described in BK, which motivate the conjecture made there, are an 
application of general techniques associated with the theory of singularity- dominated strong fluc- 
tuations. These techniques have been applied previously to analyze twinkling starlight van 
Hove-type singularities P), and the influence of classical periodic orbit bifurcations on quantum 
energy level (|| and wavefunction [|l3| statistics. In all of these applications the results correspond to 
power-law asymptotics of the moments of fluctuating quantities as the relevant parameter vanishes, 
with exponents that emerge from a competition between different singular contributions. It was 
shown by Hannay J§|, ^) for the the moments of the intensity fluctuations beyond a one-dimensional 
refracting screen that exactly when one kind of singularity overtakes another in the competition 
there is an additional logarithmic factor. Hannay also obtained the constants multiplying the 
various asymptotic contributions in this case. Importantly, in none of the applications studied pre- 
viously has it been possible to prove non-trivial predictions of the theory of singularity-dominated 
strong fluctuations by an asymptotic analysis that could be made rigorous. 

In the example we study here, the singularity competition considered in BK is between clusters 
of nearly degenerate eigenvalues. Clusters involving p eigenvalues give rise to a contribution to the 
ensemble average of \Z^{^)\~ k that diverges like e~ Vp ( k ' as e — -> 0. For a given k, the dominating 
cluster-size is the one for which the exponent v p {k) is maximal. It was shown in BK that this 
produces the exponent (Q). Here, in Section ^, we show that for k an integer, when one p takes over 
from another as dominant, there is an additional logarithmic factor, as described above. In Section 
we prove this result by calculating the GOE average explicitly, in the large matrix-size limit. 
This represents the first example where nontrivial predictions of theory of singularity-dominated 
strong fluctuations have been proved. 



asymptotics is much simpler [^TJ . 




(4) 
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2 Cluster contributions 



We here re-analyze the arguments presented in BK to recover explicitly the logarithmic factor 
anticipated there in the case of negative integer moments of characteristic polynomials of random 
matrices in the GOE. 

Denoting by M p (—k, e) the contribution from clusters of p eigenvalues (we henceforth refer to 
this as the p- cluster contribution) to the GOE average of |Z7v(/i)| _fe , where e is Im/i scaled by the 
mean eigenvalue density, equations (9) and (10) of BK may be written 



rX nX T\ P ~\ TP 4.1 \Xrn ~ Xr 

11m- 1 lln-m+1 \ tAj m ^r, 

X ~ 4 J- X ~' J.x ~~" [( X \ + e 2)( x 2 + £ 2) _ _ _ (x 2 + £ 2)] fe / 2 " 



M p (-fc,e)oc / dx x I dx 2 ... I dx p — — rj^. (5) 



To be precise, the limits of integration were given as — oo and oo in BK. This distinction will 
be important when k is an integer, and not otherwise. A finite integration range is, in fact, more 
appropriate; in the case of the circular ensembles of random matrix theory because the eigenphases 
lie in a finite interval, and in the case of the Gaussian (or similar) ensembles because the potential 
effectively limits the range in which the eigenvalues lie. 
Making the change of the variables x m = eu m gives 

m,<-*, () oc r ** f" *• • ■ ■ f" rcare-H--'---i (6) 

J-X/t J-X/e J-X/e [{ul + l){ul + l)...{ul + l)Y 12 

It was demonstrated in BK that the p-cluster contribution dominates the fcth moment when p < 
k < p + 1 . It is straightforward to see that the integral in (||) converges as e — > in the range 
p < k < p + 1. It is then asymptotically consistent to replace the limits of integration by — oo and 
oo, and the results of BK hold without change. When k is an integer, the p = k integral diverges 
and so must be treated more carefully. 



Let 



Consider first the case when p = 1: 



MA7< ! = / (hi \ / dn 2 ... I ,/„,, r n ^ inL " m|Mm U i p/2 . (7) 





j-X/e 


j-X/e 


/ dm 


1 du 2 ■ ■ 


. / du. 


J-X/e J 


'-X/e 


J-X/e 



X/e 



[(«? + l)(u3 + l)...(u2 + l)] 



dui 



which clearly diverges like log \ as e — > 0. 
Consider next the case when p = 2: 



x/ < f x/ - iu!-« a | [ x ^ A r , 



h{X/e)= I dm I du 2 . 2 , , 2 , oc / dui / du 2 , 2 , 1V 2 : 1V (9) 

J-X/e J-X/e K + 1X^2 + 1) J-X/e J-X/e K + 1) K + 1) 

This can be written as two integrals, one associated with the first term in the numerator of the 
integrand and the other associated with the second term. It may be seen straightforwardly that 
again both integrals diverge like log - as e — > 0. 
In the general case 



IpiX/e) oc / dui I du 2 ... I du„ ^1 — "~ m+ — — " ; , . (10) 

-X/e J-X/e J-X/e [( U j + 1) («| + 1) . . . (uj + 1)] ' ' 
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Expanding out the numerator of the integrand, I p may be expressed as a sum of integrals, each 
coming from a term in the resulting series. It may be seen immediately that each integral diverges 
like log - as e — > 0. Thus when k is an integer the GOE average of \Z^(^)\^ k diverges like 

e -fc(fc-D/2 log I (U) 

as e — + 0. 



3 GOE negative moments 

Our purpose now is to prove the result obtained at the end of the previous section. We shall do 
this by making a careful asymptotic analysis of the exact GOE average defining the moments. 

Regularizing the characteristic polynomial Z/v(/j) = det ^/iljv — by taking Im/i > 0, one 
may represent negative half-integer powers of the determinant as a Gaussian integral: 



-n/2i _ 



f]_ dS k exp I \ ^2 s l s k - ^Tr 



fc=i 



fc=i 



k=l 



(12) 



where we have introduced real-valued N— dimensional vectors S k = (sk,i, Sfc,Jv) T for k = 
1, 2, n so that dS k = ds k<i . 

Denoting by (...) the expectation value with respect to the distribution (|l|), our goal is to 
calculate the negative integer moments 



as well as the correlation function 

^,nG"i>M2) = ([Z N (p 1 )Z N (p* 2 )}~ n/2 



(13) 
(14) 

assuming Im(/Ji) = Im(/i2) > 0. It will be convenient for us to define fi\ — /i + % + iS and 
^2 = Li — % — i5, with fi, u> and S real and 5 > 0. Note that when lu = the correlation function 
reduces to the negative integer moments of the absolute value of the characteristic polynomial, 
which are the main objects of interest here. 

We start with (|l^) . Performing the ensemble averaging in the standard way using the identity 



/ 



dHViH^in^} = exp^-^Tr 



A 2 + AA 



(15) 



gives 



i r " I ? ™ J 2 - 

Yl dS k exp { - Ml £ S^Sfe - — £ ( S fe S <) ( S FSa 



(2iri) nN / 2 J 11 r 1 2 r "^^~ K ~" 8N 

y ' J k=l I k=l k,l=l 



(16) 



Introducing an n x n real symmetric matrix Q with matrix elements Q k i — S^S;, we note that 
the integrand may be conveniently rewritten in the form 

* _ A J 2 



exp { - Ml TrQ - — TrQ' 
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This fact allows us to employ the " integration theorem" proved in Appendix A of H and to rewrite 
the integral in ( |l6| ) in terms of an integral over the positive definite matrices Q: 

K$ n = [ dQe^l^TTQ+^TQ^} det Q(N-n-l)/2 

JQ>0 



provided N > n+ 1. We have also rescaled the integration variable: Q — > 2NQ so that the overall 



constant n is given by 



where T(z) is the Euler gamma-function. 

As the last step of the procedure we choose the eigenvalues q\,...,q n and the corresponding 
eigenvectors of Q as new integration variables. This corresponds to the change of the volume 
element 

1 " 

dQ = -G n \A{q}\T\d qi d»(O n ), (18) 
n! ■ L1 - 

i=i 

where Ajq} = Yii<j {Qi~Qj) ^ s the Vandermonde determinant and d/i(O n ) stands for the normalized 
invariant measure on the orthogonal group O(n). Here 

c " = w ^n-k) (19) 

and the factor 1/n! ensures that the integration domain with respect to all variables qk can be 
taken to be < qk < oo. 

The integrand is obviously 0(n) invariant and so we obtain: 

K&CMX) = Gj» n jf >q II (dq^'^qT^) |A{g}|exp|-|X:^fe)} - ( 2 °) 
where C^ ] n = ^G n C$ n and 

A(q) = J 2 q 2 -2ifiq-\nq. (21) 

We are mainly interested here in the limit of large matrix size, where one expects the results 
to show universality. To extract the leading asymptotics as N — * oo when n is fixed we employ 
the saddle-point method, and consider Nlo as well as NS to be of the order unity when N — * oo. 
The stationary points of A(qi) are obviously given by 

2 J 2 qi -2in =0, (22) 

ft 

where i = 1, 2, n. Each of these equations has two solutions: 



1± = • ( 23 ) 

We would like to choose the spectral parameter [i to satisfy < J\/2 in accordance with the 
idea of considering the bulk of the spectrum for GOE matrices of large size. Then only for q + are 
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the real parts positive, and so only in this case do the corresponding saddle points contribute to 
the integral over the positive semiaxis q > 0. Consequently, among the 2™ possible sets of saddle 
points (q±, ■■■,q±) only the choice 

q+ = diag(g+,...,Q + ) (24) 

is relevant. 

The presence of the Vandermonde determinants makes the integrand vanish at the saddle-point 
sets and so care should be taken when calculating the leading order contribution to the integral. 
This turns out to be given by 

/oc n 
-°° k=i fei<fe 2 

(25) 

with 

The integral in ([25| ) is a particular case of the Selberg integral |l5] and can be evaluated explicitly. 
We do not give the resulting expression here, because it is not needed for our purposes. 

We note for later purposes that a formula for /C]^ (/x^) can obviously be obtained from the 
above expression by taking its complex conjugate and then replacing [i\ with [i* 2 . 

We next consider the product of the expression (|l2|) with its complex conjugate at a different 
value of the spectral parameter and average it over the GOE. From now on we use the index 
(7 = 1,2 to label the iV-component vectors So- stemming from the first/second set of integrals. 
To write the resulting expression in a compact form it is again convenient to introduce a 2n x 2n 
matrix Q with the matrix elements Q^l'^ 2 = k^<?2,i- Here fc and / take the values 1, ...,n. In 
terms of this matrix 

1 f n { i " i n t/ 2 1 

fcffinCMi)/^) = , )Nn / ]J dSi ifc dS 2 , fc exp i ^//i^Sf ifc Si, fc - ^/X2^S| , ifc S 2 , fc - — Tr (QLQL) L 

^ ' J fc=i I fe=i fc=i J 

(27) 

where L = diag(l„, — 1„). Again employing the same integration theorem as above and changing 
Q — > 2NQ we arrive at 



/K -N.7 



(HulMi)=C%> n [ dQe-^- 2i ™$+J 2 Tr(Qm)} det Q(N-2n-i)/2^ (28) 
Jq>o 



>Q>0 

provided N > In + 1, where M = diag(/iil n , — /i|l n ) and 

= ( 7V )iVn (7r) -r 1 (2„-l)/2 



n^o lr '— 



This equation differs from its analogue (|l^) in one important aspect: it is now of little use to 
introduce the eigenvalues/eigenvectors of Q as integration variables. Rather, it is natural to treat 
Ql = QL as a new matrix to integrate over. Such (non-symmetric!) matrices satisfy Q\ — LQ^L, 
have all eigenvalues real and can be diagonalized by a (pseudo-orthogonal) similarity transformation 
Ql = TaqT^ 1 , where q = diag(<?i, —92) and the nxn diagonal matrices q\, g 2 satisfy q\ > , q 2 > 0. 
Pseudo-orthogonal matrices Tq satisfy: TqLTq — L and form the group 0(n, n) (the corresponding 



G 



symmetry is conventionally called a "hyperbolic symmetry" in the random matrix literature, see 

13). 

It turns out that a more convenient way to proceed is to block-diagonalize the matrices Q^: 

- i ( A \ - - 0(n,n) 

Ql=T- 1 [ 1 - T , where Te 



-P 2 / ' 0(n) x O(n) 

and P\ t 2 are n x n real symmetric, with positive eigenvalues 171,2, respectively. The integration mea- 
sure [<IQl\ can be derived in terms of the new variables following the standard steps (see e.g. p0[) 
outlined in the Appendix of the present paper. We arrive at [dQ] — A dP\dP 2 ]J k k2 (qi^ + <72,fe 2 ) dfi(T) , 
where A — G^/[n!2™^ rl+1 ^ 2 ] and the last factor is the invariant measure on the manifold of T- 
matrices. An explicit expression for it is presented, for reference purposes, in the Appendix. 
After all these preparatory steps we arrive at the following expression: 

Kn\ = AC Nn [ I dPldP 2 I(M,Pl,P2) (29) 
JP 1 >0JP 1 >0 



X 

ki,k 2 



I [ (qim +92,fc 2 )det 



-P1P2 



(JV-2n-l)/2 



where 



I(M, A, A) = J dn(f) exp jzATr f ^ j t 1 ( Pl _^ J t J . (30) 

Employing the explicit parametrization for the matrices T given in the Appendix we can rewrite 
the above integral as 

I(M,P U P 2 ) = e lN ^T, k ^-i^I (M,P 1 ,P 2 ), (31) 

/oa n 
n#fe n i c ° sh ^2 -coshVfe.1 (32) 
"°°fe=l ki<k 2 



x / [dfi(0 L )][diJ.(0 R )] exp i iN ^ 1 M2 Trcosh^ 



6l T Pi6 l + 6 r t P 2 6 r 

I L L 

where Ol,r € 0{n), and $ is diagonal. 

In the case of GUE matrices studied in [0] a helpful trick under similar conditions was to perform 
the (unitary) group integrals explicitly by employing the famous Itzykson-Zuber-Harish-Chandra 
integration formula M| . The lack of an analogous formula for the orthogonal group forces us to 
take a slightly different route. 

It is easy to see that the value of this integral can depend only on the eigenvalue matrices q\ and 
q 2 . Let us therefore introduce the eigenvalues (and corresponding eigenvectors) of the Hermitian 
matrices P\ > and P 2 > as the integration variables. This results in the following expression: 

poo />oo 

>4 2) >i,M2) = c% i{dq hl q-r 1/2 \M4i}\ n^g 2 :r 1/2 i A te}i 

x II (9iM +q 2M )PMAiA2)e-^ E^^(*.0 (33) 

ki ,k2 
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where 



r<( 2 ) \ r<4 /^(2) .-._> I i 



= J A q z - 2ifiq - In g and A* (<?) = J V + 2i/i§ - In g . (35) 

Again, we need to perform an asymptotic analysis as N — > oo. The most interesting regime 
occurs when one keeps the difference Re(^i — ^2) = w an d the regularization 5 so small as to ensure 
TVmax (w, S) < 00, while /i = Re ^ 1 ^ 2 ^ is kept in the range < JV2. 

The stationary points of A(q) and A* (q) are now given by 

qn~ ip, = and q 2 i + = 0, (36) 

qi.i 92, i 

where i = 1, 2, n. Each of these two equations has two solutions: 



ifi±y/2J2-^ -z M ± ^2J 2 - M 2 
Qi± = ?2± = 2J2 ' ^ > 

but only for gi + ,<72+ = <7i+ are the real parts positive; that is, only then do the corresponding 
saddle points contribute to the integral over the positive semiaxis q± t i > or qi.i > 0. Consequently, 
among the 2 2 ™ possible sets of stationary points only the choice 

qi = diag(qi+, q 1+ ) , q\ = diag(g* + , q* + ) (38) 

is relevant. This is a major simplification, because for such a choice the integrand in ( |3l| ) turns 
out to be independent of the matrices Oi, O2. 

Taking care of the Vandermonde determinants when calculating the fluctuations around the 
chosen saddle points and remembering that 

qi +q{=np(p) , q 1 q* 1 =l/2J 2 (39) 

where p(/i) = -^py/2J 2 — /i 2 is the mean density of eigenvalues for GOE matrices, we observe 
that when the asymptotic expression for the correlation function under consideration is divided 
by the product of the negative moments (^) the Selberg integrals cancel out, as well as all of the 
exponential factors too. The resulting expression amounts to 

n/2 

det (/xi In — H) det (^ 1 N — H " 



IC n ( Pl ,m) = lim ; ; ' = C x F^ OE (e) , (40) 

(det(^lN-H) n/2 )(det(»* 2 l N -H) " /2 ) 

where 

F° oE (e) = e - r n n - w 

J i k =i V A fc " 1 fcl<fc2 
in which we have introduced variables Afc = cosh^ £ [1, 00), 

e = -W 7 rp(/i)(/i 1 -/i|)/2 ) (42) 

and 

2 /iV\™ 2/2 (2ttW 2 1 

V2y n! n; =1 r(i) 
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This expression is valid for all |e| < oo, i.e. as far as (/ii — /ij) = 0(1/N) and constitutes one of 
the main results of the present paper. In the subsequent analysis we concentrate on the moments 
of characteristic polynomials and thus treat e as a real parameter. 

It is instructive to compare ( f4l| ) with its counterpart for the Gaussian Unitary Ensemble (see 
fl2lj; in H the corresponding expression is implicit): 



F, 



GUEi 



(e) = e" e / [J dX k (Afc, ~ A fe2 fe^ E*= t A * . (44) 



k=l fci</C2 



The latter integral is a specific case of the Selberg integral [|l5| and can be immediately evaluated, 
yielding 

F^ UE (e) = ^l[jKj + 1)!. (45) 

Such a formula exemplifies a 'normal' dependence of the negative moments on e: namely, one can 
extract the rate of divergence as e — * by analysing the perturbative expansion of the integral as 
e — > oo. Performing the latter limit is effectively the same as considering the case when S = Im/11,2 
is left unsealed by the mean eigenvalue density (see the Introduction). In other words, it is 
equivalent to considering the limit 6 — > after taking N — ► oo. Thus for the GUE the asymptotics 
of the negative moments is the same irrespective of the order in which limits are taken, and so is 
relatively uninteresting. 

The integral ([ll|) behaves in this sense 'anomalously'. It does not belong to the class of Selberg 
integrals and apart from when n — 1 (in which case it just yields the Macdonald function Kq(c)) 
we have failed to evaluate it explicitly in a simple closed form. We therefore proceed to analyze 
the limits e — > oo and e — > separately. 

In the perturbative region e >> 1 the integral is obviously dominated by a small vicinity of 
the lower limit: Afc — I << 1. Introducing variables x k € [0, oo) such that Afc = 1 + x^/e, we 
immediately see that asymptotically the integral is again of Selberg type: 

^ G °^»1) = 2^W /°°n% II l^-^l^-** (46) 

J I 1 V ^ 7,, ^ J„ _ 



k—1 V fci</C2 

r(¥)r(¥) 



2 n/2 £ n2/2 II p (|*) e «2/2 (2tt)"/ 2 



nr(f 

J'=l 



2 



(47) 



We then see that the perturbative behaviour for GOE moments is essentially of the same type as 
that for GUE moments: 

^ (M , M)H ^ )jr -(|)""^(-^_)"" /2 (48) 

In contrast to this, in the non-perturbative region e — > the behaviour of the GUE and GOE 
moments is very different. In this limit the integral is dominated by Afc ~ eT 1 >> I and it is 
natural to introduce rescaled variables yk = eAfc, leading to 

-i />00 />OQ n T 



Vk 

fe=l yK fci<fc 2 
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Note that one cannot set the lower limit of integration with respect to the variables to be zero, 
because the corresponding integrals diverge logarithmically there. To extract the leading order 
behaviour as e — > we differentiate the function F n (e) = e™(" -1 )/ 2 F^ OE (e « 1) with respect to 
its argument, reducing it asymptotically to a Selberg-type integral 

l&(e) = -V« r^e-ia... r^ e -«» (ltt _ cMyn _ e) J] |j, fc - w |(50) 

— -7/ II^ fc x II ly*-wl = -7ll c Vm 2 ■ ^ 

J ° fc=l l<k<2<(n-l) i=0 V27 

Thus, we conclude that for all integer n > 1 

3=0 

This anomalous behaviour paramctrically agrees with that predicted by the heuristic theory of 
dominating singularities outlined in Section I2I see in particular (|Tl|) . 



1 + | 



lnl/e 

f n(n-l)/2 ' 



(52) 
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Appendix - Calculation of the Jacobian 

To evaluate the integral in ( |28| ) one needs to calculate the Jacobian generated by the variable 
transformation Q L = f^PT, with P = diag(A, -A) = P T and P 1 > 0, P 2 > being real 
symmetric n x n matrices. For the matrices T £ o(n)x'o(n) we em Pl°y the following explicit 
parametrisation in terms of a real n x n matrix t: 

VT+W -i \ 

It is convenient to follow the scheme suggested in p0~| . One starts by considering the relation 
between the matrix differentials: 

dQ = f- 1 dQT , dQ = dP + (Pdf - dffj , (53) 

where we have introduced the notation df = dTT^ 1 . To calculate the Jacobian the difference 
between dQ and dQ is immaterial and we omit the tilde henceforth. Partitioning the matrix dQ 
into four n x n sub-blocks dq pq , p, q = 1, 2, one then rewrites the above relation blockwise: 

dq n = dPi + (Pidru ~ dfjxPj) , dq 22 = -dP 2 - (p 2 dr 22 - dT 22 P 2 ) 



f={ A—) hence P = 
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dqi2 = A df 12 + df 12 P 2 , dq 21 = dq[ 2 . 



Inspecting the block structure of the corresponding Jacobian, symbolically written as J 
dct (d [qu, q 22 , qi 2 ] jd P\,P 2 ,f ^ , one may easily verify that 

n 

J = dct (d [q 12 ] jd [f 12 ]) = dct (A (g) 1„ + 1„ ® P 2 ) = J| + ?2,j) , 



where q lti and g 2 ,i are (positive) eigenvalues of the matrices A, A- Then an intermediate result 
for the measure can be schematically written as 

n 

dQ = ]J + «2,j) dAdAdet (d [fi 2 ] /d [t] ) dt, 

where, explicitly 



-d[fi 2 ] = [dT n ] [T _1 ] 12 + [dT} 12 [T-\ 2 = d [VI + « T J i+dtVl + i T i. 

To calculate the remaining determinant we employ the singular value decomposition i = 
expressing i in terms of the two real orthogonal nxn matrices Ol,r € O(n) and a real 
diagonal matrix 9 = diag(6*i, 



assuming, for uniqueness, 6\ > ... > 6„ 



L 1 cosh OOl and + t T i = coshOOft. Further introducing dv^^R 
df = Oi,d [fi 2 ] Oft 1 we find, after straightforward manipulations, 



dO 



Then \n+W = 
1 



L,R 



O 



L,R 



and 



df = d9 + sinh 9dvR cosh 9 — cosh 9dvL sinh 0. 

Next, differentiating Ol,rO^ r = 1, we observe that dvL,R must be antisymmetric, hence diu = 
and dij^i = —diij, from which it is clear that [df],; — 0{ for all i = 1, ...,n. At the same time, for 
any of the n(n — l)/2 pairs 1 < i < j < n we have, in vector notation, the relation between the 
differentials 



dfij 
dfu 



sinh 9i cosh 9j — cosh 6i sinh #j 
- cosh 9; sinh 0, 



(dv R )ij 
[dv L )ij 



(54) 



sinh 9i cosh 0j 

The Jacobian in question then reduces to a product of the determinants of the matrices entering 

,, - | cosh 20; - cosh 20j|/2. Finally, 



in the above equation, which arc simply | sinh 2 Oi — sinh 2 9j 



we introduce ipi = 26 \ , remove the relative ordering of ipi in favour of the factor 1/n! in the 
measure, and remember that [dvL.R.) gives rise to the product of invariant measures [(1[i,(Ol,r)] on 
the orthogonal group 0(n) (which we assumed to be normalized to unity). The measure [dQi] in 
the coordinates Pi, 2 ,iP,Ol,r then assumes the following form: 



[dQ\= 2 „ ( „fi"/2 n , n(gM + g2,j) II |cosh^-cosh^|dAdP 2 ^[dM(0 L )][d^(Ofl)] (55) 

i.j l<i<j<n 



where — oo < ipi < oo for i = 1, ...n. 
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To conclude, we give the explicit expression for the following combination used in the main 
text: 



Tr 

Tr 
Tr 



Ailr 



Pi 



-Pi 



cosh 0/2 — sinh 0/2 
— sinh 0/2 cosh0/2 



cosh -0/2 sinh 0/2 
sinh 0/2 cosh 0/2 



P L (^i cosh 2 0/2-/4 sinh 2 0/2^ - Tr P R (V cosh 2 0/2 - fi sinh 2 0/2^ 



= + M2)Tr (Pl - Pfl) + - 14) ( T ? P L cosh0 + TrP fl cosh0) 



(56) 



where we have introduced matrices Pl = OlP\O l 1 and Pr = OrPzOr having the same eigen- 
values qi = diag(gi i i, ...,qi, n ) and q 2 = diag(g2,i, —,Q2, n ) as the matrices Pi i2 - 
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